Splicing analysis of STAT3 tandem donor suggests non-canonical binding registers for U1 and U6 snRNAs.

Tandem donor splice sites (5'ss) are unique regions with at least two GU dinucleotides serving as splicing cleavage sites. The Δ3 tandem 5'ss are a specific subclass of 5'ss separated by 3 nucleotides which can affect protein function by inserting/deleting a single amino acid. One 5'ss is typically preferred, yet factors governing particular 5'ss choice are not fully understood. A highly conserved exon 21 of the STAT3 gene was chosen as a model to study Δ3 tandem 5'ss splicing mechanisms. Based on multiple lines of experimental evidence, endogenous U1 snRNA most likely binds only to the upstream 5'ss. However, the downstream 5'ss is used preferentially, and the splice site choice is not dependent on the exact U1 snRNA binding position. Downstream 5'ss usage was sensitive to exact nucleotide composition and dependent on the presence of downstream regulatory region. The downstream 5'ss usage could be best explained by two novel interactions with endogenous U6 snRNA. U6 snRNA enables the downstream 5'ss usage in STAT3 exon 21 by two mechanisms: (i) binding in a novel non-canonical register and (ii) establishing extended Watson-Crick base pairing with the downstream regulatory region. This study suggests that U6:5'ss interaction is more flexible than previously thought.

Genetic findings in Czech patients with limb girdle muscular dystrophy.

Limb girdle muscular dystrophies (LGMD) are a genetically heterogeneous group of muscular dystrophies. The study presents an overview of molecular characteristics of a large cohort of LGMD patients who are representative of the Czech LGMD population. We present 226 LGMD probands in which 433 mutant alleles carrying 157 different variants with a supposed pathogenic effect were identified. Fifty-four variants have been described only in the Czech LGMD population so far. LGMD R1 caplain3-related is the most frequent subtype of LGMD involving 53.1% of patients with genetically confirmed LGMD, followed by LGMD R9 FKRP-related (11.1%), and LGMD R12 anoctamin5-related (7.1%). If we consider identified variants, then all but five were small-scale variants. One large gene deletion was identified in the LAMA2 gene and two deletions in each of CAPN3 and SGCG. We performed comparison our result with other published studies. The results obtained in the Czech LGMD population clearly differ from the outcome of other LGMD populations in two aspects-we have a more significant proportion of patients with LGMD R1 calpain3-related and a smaller proportion of LGMD R2 dysferlin-related.

Myeloid lineage cells evince distinct steady-state level of certain gene groups in dependence on hereditary angioedema severity.

Hereditary angioedema (HAE) is a rare genetic disorder with variable expressivity even in carriers of the same underlying genetic defect, suggesting other genetic and epigenetic factors participate in modifying HAE severity. Recent knowledge indicates the role of immune cells in several aspects of HAE pathogenesis, which makes monocytes and macrophages candidates to mediate these effects. Here we combined a search for HAE phenotype modifying gene variants with the characterization of selected genes' mRNA levels in monocyte and macrophages in a symptom-free period. While no such gene variant was found to be associated with a more severe or milder disease, patients revealed a higher number of dysregulated genes and their expression profile was significantly altered, which was typically manifested by changes in individual gene expression or by strengthened or weakened relations in mutually co-expressed gene groups, depending on HAE severity. SERPING1 showed decreased expression in HAE-C1INH patients, but this effect was significant only in patients carrying mutations supposedly activating nonsense-mediated decay. Pro-inflammatory CXC chemokine superfamily members CXCL8, 10 and 11 were downregulated, while other genes such as FCGR1A, or long non-coding RNA NEAT1 were upregulated in patients. Co-expression within some gene groups (such as an NF-kappaB function related group) was strengthened in patients with a severe and/or mild course compared to controls. All these findings show that transcript levels in myeloid cells achieve different activation or depression levels in HAE-C1INH patients than in healthy controls and/or based on disease severity and could participate in determining the HAE phenotype.

RNA kink-turns are highly anisotropic with respect to lateral displacement of the flanking stems.

Kink-turns are highly bent internal loop motifs commonly found in the ribosome and other RNA complexes. They frequently act as binding sites for proteins and mediate tertiary interactions in larger RNA structures. Kink-turns have been a topic of intense research, but their elastic properties in the folded state are still poorly understood. Here we use extensive all-atom molecular dynamics simulations to parameterize a model of kink-turn in which the two flanking helical stems are represented by effective rigid bodies. Time series of the full set of six interhelical coordinates enable us to extract minimum energy shapes and harmonic stiffness constants for kink-turns from different RNA functional classes. The analysis suggests that kink-turns exhibit isotropic bending stiffness but are highly anisotropic with respect to lateral displacement of the stems. The most flexible lateral displacement mode is perpendicular to the plane of the static bend. These results may help understand the structural adaptation and mechanical signal transmission by kink-turns in complex natural and artificial RNA structures.

Next-generation sequencing in children with epilepsy: The importance of precise genotype-phenotype correlation.

AIM: The primary goal was to determine the yield of next-generation sequencing (NGS) epilepsy gene panels used for epilepsy etiology diagnosing using a multidisciplinary approach and to demonstrate the importance of genotype-phenotype correlations. The secondary goal was to evaluate the application of precision medicine in selected patients. METHODS: This single-center retrospective study included a total of 175 patients (95 males and 80 females) aged 0-19 years. They were examined between 2015 and 2020 using an NGS epilepsy gene panel (270 genes). A bioinformatic analysis was performed including copy number variation identification. Thorough genotype-phenotype correlation was performed. RESULTS: Out of 175 patients, described pathogenic variants or novel variants with clear pathogenic impact were identified in 30 patients (17.14%). Genotype-phenotype correlations and parental DNA analysis were performed, and genetic diagnosis was confirmed on the basis of the results in another 16 out of 175 patients (9.14%). The diagnostic yield of our study increased from 30 to 46 patients (by 53.33%) by the precise genotype-phenotype correlation. INTERPRETATION: We emphasize a complex genotype-phenotype correlation and a multidisciplinary approach in evaluating the results of the NGS epilepsy gene panel, which enables the most accurate genetic diagnosis and correct interpretation of results.

A GP1BA Variant in a Czech Family with Monoallelic Bernard-Soulier Syndrome.

Bernard-Soulier syndrome (BSS) is a rare inherited disorder characterized by unusually large platelets, low platelet count, and prolonged bleeding time. BSS is usually inherited in an autosomal recessive (AR) mode of inheritance due to a deficiency of the GPIb-IX-V complex also known as the von Willebrand factor (VWF) receptor. We investigated a family with macrothrombocytopenia, a mild bleeding tendency, slightly lowered platelet aggregation tests, and suspected autosomal dominant (AD) inheritance. We have detected a heterozygous GP1BA likely pathogenic variant, causing monoallelic BSS. A germline GP1BA gene variant (NM_000173:c.98G > A:p.C33Y), segregating with the macrothrombocytopenia, was detected by whole-exome sequencing. In silico analysis of the protein structure of the novel GPIbα variant revealed a potential structural defect, which could impact proper protein folding and subsequent binding to VWF. Flow cytometry, immunoblot, and electron microscopy demonstrated further differences between p.C33Y GP1BA carriers and healthy controls. Here, we provide a detailed insight into its clinical presentation and phenotype. Moreover, the here described case first presents an mBSS patient with two previous ischemic strokes.

Phylogenetic Reassessment, Taxonomy, and Biogeography of Codinaea and Similar Fungi.

The genus Codinaea is a phialidic, dematiaceous hyphomycete known for its intriguing morphology and turbulent taxonomic history. This polyphasic study represents a new, comprehensive view on the taxonomy, systematics, and biogeography of Codinaea and its relatives. Phylogenetic analyses of three nuclear loci confirmed that Codinaea is polyphyletic. The generic concept was emended; it includes four morphotypes that contribute to its morphological complexity. Ancestral inference showed that the evolution of some traits is correlated and that these traits previously used to delimit taxa at the generic level occur in species that were shown to be congeneric. Five lineages of Codinaea-like fungi were recognized and introduced as new genera: Codinaeella, Nimesporella, Stilbochaeta, Tainosphaeriella, and Xyladelphia. Dual DNA barcoding facilitated identification at the species level. Codinaea and its segregates thrive on decaying plants, rarely occurring as endophytes or plant pathogens. Environmental ITS sequences indicate that they are common in bulk soil. The geographic distribution found using GlobalFungi database was consistent with known data. Most species are distributed in either the Holarctic realm or tropical geographic regions. The ancestral climatic zone was temperate, followed by transitions to the tropics; these fungi evolved primarily in Eurasia and Americas, with subsequent transitions to Africa and Australasia.

Functional analysis of germline ETV6 W380R mutation causing inherited thrombocytopenia and secondary acute lymphoblastic leukemia or essential thrombocythemia.

Germline mutations in ETV6 gene cause inherited thrombocytopenia with leukemia predisposition. Here, we report on functional validation of ETV6 W380R mutation segregating with thrombocytopenia in a family where two family members also suffered from acute lymphoblastic leukemia (ALL) or essential thrombocythemia (ET). In-silico analysis predicted impaired DNA binding due to W380R mutation. Functional analysis showed that this mutation prevents the ETV6 protein from localizing into the cell nucleus and impairs the transcriptional repression activity of ETV6. Based on the germline ETV6 mutation, ET probably started with somatic JAK2 V617F mutation, whereas ALL could be caused by diverse mechanisms: high-hyperdiploidity; somatic deletion of exon 1 IKZF1 gene; or somatic mutations of other genes found by exome sequencing of the ALL sample taken at the diagnosis.

Importance of base-pair opening for mismatch recognition.

Mismatch repair is a highly conserved cellular pathway responsible for repairing mismatched dsDNA. Errors are detected by the MutS enzyme, which most likely senses altered mechanical property of damaged dsDNA rather than a specific molecular pattern. While the curved shape of dsDNA in crystallographic MutS/DNA structures suggests the role of DNA bending, the theoretical support is not fully convincing. Here, we present a computational study focused on a base-pair opening into the minor groove, a specific base-pair motion observed upon interaction with MutS. Propensities for the opening were evaluated in terms of two base-pair parameters: Opening and Shear. We tested all possible base pairs in anti/anti, anti/syn and syn/anti orientations and found clear discrimination between mismatches and canonical base-pairs only for the opening into the minor groove. Besides, the discrimination gap was also confirmed in hotspot and coldspot sequences, indicating that the opening could play a more significant role in the mismatch recognition than previously recognized. Our findings can be helpful for a better understanding of sequence-dependent mutability. Further, detailed structural characterization of mismatches can serve for designing anti-cancer drugs targeting mismatched base pairs.

Low Density Lipoprotein Receptor Variants in the Beta-Propeller Subdomain and Their Functional Impact.

Background: Pathogenic variants in the low density lipoprotein receptor gene are associated with familial hypercholesterolemia. Some of these variants can result in incorrect folding of the LDLR protein, which is then accumulated inside the cell and cannot fulfill its function to internalize LDL particles. We analyzed the functional impact of 10 LDLR variants localized in the beta-propeller of epidermal growth factor precursor homology domain. The experimental part of the work was complemented by a structural analysis on the basis of 3D LDLR protein structure. Methods: T-Rex Chinese hamster ovary cells transfected with the human LDLR gene were used for live cell imaging microscopy, flow cytometry, and qRT-PCR analysis. Results: Our results showed that the analyzed LDLR protein variants can be divided into three groups. (1) The variants buried inside the 3D protein structure expressing proteins accumulated in the endoplasmic reticulum (ER) with no or reduced plasma membrane localization and LDL particle internalization, and associated with an increased gene expression of ER-resident chaperones. (2) The variants localized on the surface of 3D protein structure with slightly reduced LDLR plasma membrane localization and LDL particle internalization, and associated with no increased mRNA level of ER-resident chaperones. (3) The variants localized on the surface of the 3D protein structure but expressing proteins with cell responses similar to the group 1. Conclusion: All analyzed LDLR variants have been evaluated as pathogenic but with different effects on protein localization and function, and expression of genes associated with ER stress.

Bending of DNA duplexes with mutation motifs.

Mutations can be induced by environmental factors but also arise spontaneously during DNA replication or due to deamination of methylated cytosines at CpG dinucleotides. Sites where mutations occur with higher frequency than would be expected by chance are termed hotspots while sites that contain mutations rarely are termed coldspots. Mutations are permanently scanned and repaired by repair systems. Among them, the mismatch repair targets base pair mismatches, which are discriminated from canonical base pairs by probing altered elasticity of DNA. Using biased molecular dynamics simulations, we investigated the elasticity of coldspots and hotspots motifs detected in human genes associated with inherited disorders, and also of motifs with Czech population hotspots and de novo mutations. Main attention was paid to mutations leading to G/T and A+/C pairs. We observed that hotspots without CpG/CpHpG sequences are less flexible than coldspots, which indicates that flexible sequences are more effectively repaired. In contrary, hotspots with CpG/CpHpG sequences exhibited increased flexibility as coldspots. Their mutability is more likely related to spontaneous deamination of methylated cytosines leading to C > T mutations, which are primarily targeted by base excision repair. We corroborated conclusions based on computer simulations by measuring melting curves of hotspots and coldspots containing G/T mismatch.

Structural and Functional Impact of Seven Missense Variants of Phenylalanine Hydroxylase.

The molecular genetics of well-characterized inherited diseases, such as phenylketonuria (PKU) and hyperphenylalaninemia (HPA) predominantly caused by mutations in the phenylalanine hydroxylase (PAH) gene, is often complicated by the identification of many novel variants, often with no obvious impact on the associated disorder. To date, more than 1100 PAH variants have been identified of which a substantial portion have unknown clinical significance. In this work, we study the functionality of seven yet uncharacterized PAH missense variants p.Asn167Tyr, p.Thr200Asn, p.Asp229Gly, p.Gly239Ala, p.Phe263Ser, p.Ala342Pro, and p.Ile406Met first identified in the Czech PKU/HPA patients. From all tested variants, three of them, namely p.Asn167Tyr, p.Thr200Asn, and p.Ile406Met, exerted residual enzymatic activity in vitro similar to wild type (WT) PAH, however, when expressed in HepG2 cells, their protein level reached a maximum of 72.1% ± 4.9%, 11.2% ± 4.2%, and 36.6% ± 7.3% compared to WT PAH, respectively. Remaining variants were null with no enzyme activity and decreased protein levels in HepG2 cells. The chaperone-like effect of applied BH4 precursor increased protein level significantly for p.Asn167Tyr, p.Asp229Gly, p.Ala342Pro, and p.Ile406Met. Taken together, our results of functional characterization in combination with in silico prediction suggest that while p.Asn167Tyr, p.Thr200Asn, and p.Ile406Met PAH variants have a mild impact on the protein, p.Asp229Gly, p.Gly239Ala, p.Phe263Ser, and p.Ala342Pro severely affect protein structure and function.

High-throughput analysis revealed mutations' diverging effects on SMN1 exon 7 splicing.

Splicing-affecting mutations can disrupt gene function by altering the transcript assembly. To ascertain splicing dysregulation principles, we modified a minigene assay for the parallel high-throughput evaluation of different mutations by next-generation sequencing. In our model system, all exonic and six intronic positions of the SMN1 gene's exon 7 were mutated to all possible nucleotide variants, which amounted to 180 unique single-nucleotide mutants and 470 double mutants. The mutations resulted in a wide range of splicing aberrations. Exonic splicing-affecting mutations resulted either in substantial exon skipping, supposedly driven by predicted exonic splicing silencer or cryptic donor splice site (5'ss) and de novo 5'ss strengthening and use. On the other hand, a single disruption of exonic splicing enhancer was not sufficient to cause major exon skipping, suggesting these elements can be substituted during exon recognition. While disrupting the acceptor splice site led only to exon skipping, some 5'ss mutations potentiated the use of three different cryptic 5'ss. Generally, single mutations supporting cryptic 5'ss use displayed better pre-mRNA/U1 snRNA duplex stability and increased splicing regulatory element strength across the original 5'ss. Analyzing double mutants supported the predominating splicing regulatory elements' effect, but U1 snRNA binding could contribute to the global balance of splicing isoforms. Based on these findings, we suggest that creating a new splicing enhancer across the mutated 5'ss can be one of the main factors driving cryptic 5'ss use.

Inherited ichthyoses: molecular causes of the disease in Czech patients.

Inherited ichthyoses belong to a large and heterogeneous group of mendelian disorders of cornification, and can be distinguished by the quality and distribution of scaling and hyperkeratosis, by other dermatologic and extracutaneous involvement, and by inheritance. We present the genetic analysis results of probands with X-linked ichthyosis, autosomal recessive congenital ichthyosis, keratinopathic ichthyosis, and a patient with Netherton syndrome. Genetic diagnostics was complemented by in silico missense variant analysis based on 3D protein structures and commonly used prediction programs to compare the yields of these two approaches to each other. This analysis revealed various structural defects in proteins coded by mutated genes while no defects were associated with known polymorphisms. Two patients with pathogenic variants in the ABCA12 gene have a premature termination codon mutation on one allele and a silent variant on the second. The silent variants c.69G > A and c.4977G > A are localised in the last nucleotide of exon 1 and exon 32, respectively, and probably affect mRNA splicing. The phenotype of both patients is very severe, including a picture harlequin foetus after birth; later (at 3 and 6 years of age, respectively) ectropin, eclabion, generalised large polygonal scaling and erythema.

A novel germline mutation of the SFTPA1 gene in familial interstitial pneumonia.

Different genes related to alveolar stability have been associated with familial interstitial pneumonia (FIP). Here, we report a novel, rare SFTPA1 variant in a family with idiopathic interstitial pneumonia (IIP). We performed whole-exome sequencing on germline DNA samples from four members of one family; three of them showed signs of pulmonary fibrosis (idiopathic interstitial pneumonia) with autosomal-dominant inheritance. A heterozygous single nucleotide variant c.532 G > A in the SFTPA1 gene has been identified. This variant encodes the substitution p.(Val178Met), localized within the carbohydrate recognition domain of surfactant protein A and segregates with the genes causing idiopathic interstitial pneumonia. This rare variant has not been previously reported. We also analyzed the detected sequence variant in the protein structure in silico. The replacement of valine by the larger methionine inside the protein may cause a disruption in the protein structure. The c.532 G > A variant was further validated using Sanger sequencing of the amplicons, confirming the diagnosis in all symptomatic family members. Moreover, this variant was also found by Sanger sequencing in one other symptomatic family member and one young asymptomatic family member. The autosomal-dominant inheritance, the family history of IIP, and the evidence of a mutation occurring in part of the SFTPA1 gene all suggest a novel variant that causes FIP.

Activation-induced deaminase and its splice variants associate with trisomy 12 in chronic lymphocytic leukemia.

Activation-induced cytidine deaminase (AID) is a mutator enzyme essential for somatic hypermutation (SHM) and class switch recombination (CSR) during effective adaptive immune responses. Its aberrant expression and activity have been detected in lymphomas, leukemias, and solid tumors. In chronic lymphocytic leukemia (CLL) increased expression of alternatively spliced AID variants has been documented. We used real-time RT-PCR to quantify the expression of AID and its alternatively spliced transcripts (AIDΔE4a, AIDΔE4, AIDivs3, and AIDΔE3E4) in 149 CLL patients and correlated this expression to prognostic markers including recurrent chromosomal aberrations, the presence of complex karyotype, mutation status of the immunoglobulin heavy chain variable gene, and recurrent mutations. We report a previously unappreciated association between higher AID transcript levels and trisomy of chromosome 12. Functional analysis of AID splice variants revealed loss of their activity with respect to SHM, CSR, and induction of double-strand DNA breaks. In silico modeling provided insight into the molecular interactions and structural dynamics of wild-type AID and a shortened AID variant closely resembling AIDΔE4, confirming its loss-of-function phenotype.

A novel germline mutation in GP1BA gene N-terminal domain in monoallelic Bernard-Soulier syndrome.

Mutations in the GP1BA gene have been associated with platelet-type von Willebrand disease and Bernard-Soulier syndrome. Here, we report a novel GP1BA mutation in a family with autosomal dominant macrothrombocytopenia and mild bleeding. We performed analyses of seven family members. Using whole-exome sequencing of germline DNA samples, we identified a heterozygous single-nucleotide change in GP1BA (exone2:c.176T>G), encoding a p.Leu59Arg substitution in the N-terminal domain, segregating with macrothrombocytopenia. This variant has not been previously reported. We also analysed the structure of the detected sequence variant in silico. In particular, we used the crystal structure of the human platelet receptor GP Ibα N-terminal domain. Replacement of aliphatic amino-acid Leu 59 with charged, polar and larger arginine probably disrupts the protein structure. An autosomal dominant mode of inheritance, a family history of mild bleeding episodes, aggregation pattern in affected individuals together with evidence of mutation occurring in part of the GP1BA gene encoding the leucine-rich repeat region suggest a novel variant causing monoallelic Bernard-Soulier syndrome.

C-terminal RUNX1 mutation in familial platelet disorder with predisposition to myeloid malignancies.

Here we report a C-terminal RUNX1 mutation in a family with platelet disorder and predisposition to myeloid malignancies. We identified the mutation c.866delG:p.Gly289Aspfs*22 (NM_001754) (RUNX1 b-isoform NM_001001890; c.785delG:p.Gly262Aspfs*22) using exome sequencing of samples obtained from eight members of a single family. The mutation found in our pedigree is within exon eight and the transactivation domain of RUNX1. One of the affected individuals developed myelodysplastic syndrome (MDS), which progressed to acute myelogenous leukemia (AML). A search for the second hit which led to the development of MDS and later AML in this individual revealed the PHF6 gene variant (exon9:c.872G > A:p.G291E; NM_001015877), BCORL1 (exon3:c.1111A > C:p.T371P; NM_001184772) and BCOR gene variant (exon4:c.2076dupT:p.P693fs; NM_001123383), which appear to be very likely second hits participating in the progression to myeloid malignancy.

Role of Inosine⁻Uracil Base Pairs in the Canonical RNA Duplexes.

Adenosine to inosine (A⁻I) editing is the most common modification of double-stranded RNA (dsRNA). This change is mediated by adenosine deaminases acting on RNA (ADARs) enzymes with a preference of U>A>C>G for 5′ neighbor and G>C=A>U or G>C>U=A for 3′ neighbor. A⁻I editing occurs most frequently in the non-coding regions containing repetitive elements such as ALUs. It leads to disruption of RNA duplex structure, which prevents induction of innate immune response. We employed standard and biased molecular dynamics (MD) simulations to analyze the behavior of RNA duplexes with single and tandem inosine⁻uracil (I⁻U) base pairs in different sequence context. Our analysis showed that the I⁻U pairs induce changes in base pair and base pair step parameters and have different dynamics when compared with standard canonical base pairs. In particular, the first I⁻U pair from tandem I⁻U/I⁻U systems exhibited increased dynamics depending on its neighboring 5′ base. We discovered that UII sequence, which is frequently edited, has lower flexibility compared with other sequences (AII, GII, CII), hence it only modestly disrupts dsRNA. This might indicate that the UAA motifs in ALUs do not have to be sufficiently effective in preventing immune signaling.